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Exemples  d’utilisation  des  techniques  d’optimisation 
en  calcul  de  structures  de  reacteurs 


Christian  Paleczny  et  David  Monteiro-Fernandes 

Snecma  Moteurs 
Centre  de  Villaroche 
77550  Moissy  Cramayel 
France 


Introduction 

Le  contexte  general  de  la  competition  sur  le  marche  militaire  de  I’aeronautique  oblige  les  differents  acteurs  a 
reduire  leurs  couts  et  a augmenter  les  performances  de  leurs  produits.  Cela  implique  un  effort  dans  tons  les 
domaines  qui  doit  se  traduire  notamment  par  la  reduction  des  cycles  de  developpement,  la  diminution  du 
nombre  d’essais  et  la  reduction  du  cout  de  fabrication. 

Dans  cette  optique,  les  outils  de  simulation  sont  devenus  indispensables  pour  predire  le  comportement  des 
pieces.  Cela  a permis  de  diminuer  le  nombre  d’essais  de  maniere  significative,  tout  en  concevant  des  pieces  de 
plus  en  plus  performantes  et  fiables. 

L’ evolution  importante  des  capacites  des  ordinateurs  permet  desormais  d’envisager  des  analyses  de  plus  en 
plus  complexes  et  de  plus  en  plus  nombreuses  en  phase  de  conception.  Cela  permet  egalement  de  diminuer 
encore  plus  les  risques  technologiques. 

Le  dimensionnement  mecanique  des  pieces  et  des  structures  pese  pour  une  lourde  part  dans  ces  couts  de 
developpements.  C’est  le  dimensionnement  mecanique  qui  garantit  I’integrite  physique  du  moteur. 

A ce  titre,  les  concepteurs  sont  tenus  a une  grande  rigueur  dans  le  travail  effectue  et  au  respect  des  procedures 
qualite.  L’ optimisation,  c’est-a-dire  I’atteinte  d’un  objectif  donne  (un  niveau  de  contrainte,  une  frequence) 
tout  en  respectant  un  certain  nombre  de  criteres  (variation  de  forme,  duree  de  vie,  ...)  est  un  domaine  qui 
prend  une  importance  croissante  dans  tons  les  bureaux  d’etudes  aeronautiques. 

En  optimisation  de  structures,  nous  distinguons  differentes  categories  : 

□ optimisation  a base  de  connaissances  (avant-projets) 

□ optimisation  topologique  (forme  non  definie  par  une  CAO) 

□ optimisation  de  dimensionnement  (epaisseurs  de  coques,  caracteristiques  materiau) 

□ optimisation  geometrique  (architecture  figee) 

A la  difference  du  secteur  automobile  et  des  avionneurs,  la  plupart  des  composants  des  reacteurs  necessitent 
une  optimisation  geometrique  des  le  debut  de  la  phase  de  conception.  L’ optimisation  geometrique  est  une 
categoric  assez  complexe  car  sa  mise  en  oeuvre  necessite,  dans  un  bureau  d’etudes,  le  couplage  automatique 
d’un  logiciel  de  CAO  parametree,  d’un  logiciel  de  maillages,  d’un  solveur  element  finis  et  d’un  logiciel 
d’optimisation.  Certains  logiciels  integrent  toutes  ces  fonctionnalites  de  faqon  « transparente  » et  ont  done  des 
performances  ties  interessantes.  Neanmoins,  ils  ont  des  limites  que  les  modelisations  des  composants  des 
reacteurs  atteignent  de  faqon  ties  frequente. 

C’est  pourquoi,  ces  dernieres  annees,  se  sont  developpes  des  logiciels  d’optimisation  ouverts  permettant  de 
gerer  tout  type  de  modelisation.  Ces  logiciels  apportent  des  possibilites  de  couplage  avec  des  logiciels  de 
CAO,  des  solveurs  elements  finis  ou  des  logiciels  developpes  en  interne. 

Apres  un  bref  chapitre  detaillant  les  principales  etapes  permettant  de  batir  un  probleme  d’optimisation,  cet 
article  presentera  quelques  exemples  d’utilisation  de  ces  logiciels  dans  un  contexte  de  bureau  d’etudes 
mecaniques. 


Communication  presentee  tors  du  symposium  RTOAVT  sur  « La  reduction  des  couts  et  des  delais  d’ acquisition 
des  vehicules  militaires  par  la  modelisation  avancee  et  la  simulation  de  produit  virtuel  », 
organise  a Paris,  en  France,  du  22  au  25  avril  2002,  et  editee  dans  RTO-MP-089. 
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Processus  d’optimisation 

Le  dimensionnement  classique  dans  un  bureau  d’etudes  repose  encore  souvent  sur  une  approche  essais- 
erreurs.  C’est  a dire,  pour  realiser  un  composant,  le  concepteur  doit: 

□ realiser  un  modele  CAO  initial  a partir  d’un  cahier  des  charges,  de  son  experience  et  de 
I’environnement  general  oil  se  situe  sa  piece 

□ realiser  ensuite  des  modelisations  (statique,  dynamique,  etc)  permettant  de  verifier  les  criteres  de 
dimensionnement 

□ iterer  en  modifiant  sa  conception  de  maniere  a respecter  son  cahier  des  charges 

Le  principe  de  base  des  logiciels  d’optimisation  est  de  pouvoir  reproduce  de  faqon  automatique  ce  qu’un 
concepteur  realisait  auparavant  manuellement,  en  y ajoutant  de  nombreuses  plus-values: 

□ possibilite  de  balayer  un  espace  de  conception  plus  large 

□ calculs  automatiques 

□ possibilites  de  realiser  des  plans  d’experiences  et  de  creer  ainsi  des  fonctions  d’ approximation 

□ atteinte  d’ optimum  grace  a des  algorithmes  de  plus  en  plus  performants 

Pour  le  concepteur,  les  etapes-clefs  en  optimisation  de  forme  sont  done: 

□ la  realisation  du  modele  de  CAO  parametree 

□ la  determination  des  analyses  a effectuer  (statique,  dynamique,  non  lineaire,  etc) 

□ la  realisation  d’un  modele  parametre  elements  finis  et  associe  a la  geometric 

□ la  definition  des  criteres  de  dimensionnement  associes  a la  pieces  (contraintes  maxi,  plage  de 
frequence  interdite,  masse,  deplacement  maxi,  etc) 


Figure  3 - Resultats 


Le  logiciel  d’optimisation  va  etre  charge  de  reconnaitre  et  de  gerer  ces  modeles  (Elements  Finis,  CAO,  etc), 
de  permettre  a I’utilisateur  d’ identifier  les  parametres  sensibles,  de  piloter  automatiquement  le  processus 
complet  de  calcul  et  de  traiter  les  resultats  en  fonction  des  parametres. 


Le  plan  d’experiences  se  contentera  de  permettre  la  creation  de  fonction  d’ approximation  des  fonctions  (mais 
possibilites  d’erreurs  importante  entre  la  fonction  estimee  et  les  points  reels),  tandis  que  1’ optimisation  directe 
utilisera  les  valeurs  reelles  des  fonctions  et  en  general  les  sensibilites  (influence  des  parametres  sur  les 
fonctions  = derivee  des  fonctions  par  rapport  aux  parametres)  pour  atteindre  des  valeurs  optimales. 

D’autres  utilisations  sont  possibles  : analyses  probabilistes  ou  de  Monte-Carlo,  recalage  de  modeles  (calculs- 
essais),  etc. 
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Sans  entrer  dans  les  details  techniques  (notion  de  perturbateur  de  maillage,  gestion  du  calcul  d’erreur,  etc), 
dans  le  cadre  de  I’optimisation  de  structures,  le  logiciel  d’ optimisation  va  effectuer  en  batch  le  processus 
iteratif  suivant : 
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Optimisation  d’un  disque  a I’eclatement 

Le  premier  exemple  illustre  un  cas  que  les  logiciels  d’ optimisation  de  base  ne  savent  pas  traiter,  car  il 
necessite  d’utiliser  un  code  de  mecanique  non  lineaire  et  porte  sur  un  critere  (rupture  de  piece)  rarement 
implemente.  L’exemple  a servi  dans  le  passe  a valider  I’utilisation  de  plan  d’experiences  et  d‘etudes 
d’ optimisation  au  sein  du  bureau  d’etudes. 

La  piece  concernee  est  un  disque  de  reacteur.  Le  role  essentiel  des  disques  est  la  retention  des  aubes  de 
compresseur  ou  de  turbine  et  le  maintien  des  jeux  en  bout  d’aube.  Ce  sont  des  pieces  tournantes  « a 
defaillance  catastrophique  » pour  Lavion.  Le  dimensionnement  de  ces  disques  requiert  done  le  respect  des 
regies  imposees  par  les  autorites  certifiantes  et  le  respect  d’un  ensemble  de  criteres  issu  du  savoir-faire  de 
Sneema  Moteurs. 


Figure  4 - modele  de  disque  Figure  5 - determination  de  la  vitesse  d’eclatement 


L’ experience  a montre  que  les  disques  eclatent  des  que  Ton  depasse  une  vitesse  de  rotation  critique,  a cause 
de  I’augmentation  rapide  des  charges  centrifuges.  Cette  vitesse  d’eclatement  pent  etre  obtenue  par  essai,  par 
formules  analytiques  ou  estimee  par  des  modeles  elements  finis.  A litre  d’information,  le  modele  numerique 
que  Ton  realise  utilise  des  lois  de  comportement  materiaux  non  lineaires  ainsi  que  des  hypotheses  de  grands 
deplacements. 

L’ etude  d’ optimisation  consiste  a minimiser  la  masse  du  disque  tout  en  respectant  une  vitesse  de  rupture 
superieure  a une  limite  imposee  (ex  : 1600  rad/s).  Pour  cette  etude,  nous  avons  modelise  un  disque  simple  et 
decide  de  selectionner  trois  parametres  qui  influent  sur  la  rupture  : I’epaisseur  de  la  toile  du  disque  El,  la 
hauteur  L3  et  la  largeur  L2  d’alesage.  Le  modele  realise  est  un  modele  2D.  L’aube  est  modelisee  par  un 
ensemble  de  masses  ponctuelles  situees  au  niveau  de  la  jante. 


Figure  6 - parametrage  et  maillage  du  disque 
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Le  probleme  d’ optimisation  que  Ton  se  propose  d’etudier  est  done  : 

□ minimiser  la  masse  du  disque, 

□ verifier  que  la  vitesse  d’eclatement  de  la  piece  est  superieure  a la  vitesse  limite. 


Afin  d’explorer  I’espace  de  conception  et  d’etudier  le  comportement  des  fonctions,  une  etude  parametrique  a 
ete  menee.  Cette  etude  permet  de  determiner  F evolution  de  la  vitesse  d’eclatement  et  de  la  masse  en  fonction 
des  parametres.  Cette  etude  permet  aussi  de  trouver  une  position  approchee  de  I’optimum  et  de  valider 
I’optimum  trouve  par  la  methode  directe. 


Un  exemple  d’application  numerique  est  donne  ci- 
apres.  L’espace  de  conception  est  defini  par  les 
plages  de  variation  des  parametres  suivantes: 

2 *E1 : 4.4  a 6 mm  par  pas  de  0.4  mm  (5  valeurs), 
2*L2  : 30  a 40  mm  par  pas  de  2.5  mm  f5  valeurs), 

L3  : 15  a 20  mm  par  pas  de  1.25  mm  (5  valeurs). 


Toutes  les  combinaisons  des  parametres  ont  ete  etudiees,  soit  125  calculs  non  lineaires  dans  un  domaine  de 
conception  coherent  avec  Fobjectif.  Une  fois  ces  calculs  obtenus,  nous  avons  cree  une  surface  de  reponse, 
nous  pouvons  done  avoir  une  representation  de  Fevolution  des  criteres  dimensionnants  en  fonction  des 
parametres  geometriques. 

La  meme  etude  a ete  ensuite  lancee  en  vue  d’une  optimisation  directe  en  utilisant  un  algorithme  classique 
(linearisation  convexe).  Le  logiciel  a converge  en  une  dizaine  d’ iteration  vers  un  point  optimum  tres  proche  de 
ce  que  Fon  trouve  en  analysant  les  resultats  du  plan  d’experiences. 


Figure  8 - Surface  de  reponse  et  optimum  direct 
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Optimisation  de  brides  festonnees 

Un  reacteur  est  constitue  de  nombreuses  brides  de  liaison.  Elies  sont  presentes  aussi  bien  sur  des  pieces  en 
rotation  que  sur  des  carters  structuraux. 

L’ operation  de  festonnage  des  brides  permet  de  gagner  de  la  masse,  ainsi  que  de  diminuer  les  niveaux  de 
contrainte  au  fond  des  trous  de  vis  (contrainte  generee  par  le  flux  des  contraintes  tangentielles).  Le  bon 
dimensionnement  de  la  geometrie  des  festons  permet  d’absorber  une  partie  de  ces  contraintes. 

L’ etude  que  nous  presentons  a permis  : 

□ de  valider  la  chaine  complete  d’optimisation  integrant  le  logiciel  de  CAO  (CATIA),  le  logiciel  de 
maillage  (PATRAN)  et  le  solveur  elements  finis  (SAMCEF) 

□ de  mieux  apprehender  la  problematique  « qualite  de  maillage  » et  « optimisation  »,  lorsque  Ton  ne 
dispose  pas  de  techniques  de  « perturbateur  de  maillage  » (re-localisation  automatique  des  noeuds  du 
maillage  en  fonction  des  variations  de  geometrie),  ni  de  calculs  de  sensibilites  semi-analytiques. 

Le  modele  de  bride  parametree  est  decrit  dans  les  figures  suivantes.  En  utilisant  les  symetries,  on  ne  se 
preoccupe  que  d’une  demi-periode  circulaire  de  cette  bride. 


Get  exemple  permet  de  montrer  les  differences  de  resultats  a travers  differentes  strategies  de  maillage  et  done 
d’illustrer  f impact  sur  une  eventuelle  optimisation  de  la  bride.  En  effet,  I’optimisation  necessite  le  calcul  des 
sensibilites  qui,  dans  ce  cas,  sont  calculees  par  differences  fmies. 


Trois  strategies  sont  etudiees  : 


□ un  maillage  libre  (figure  11), 

□ un  maillage  regie  (figure  12), 

□ un  maillage  hybride  (figure  13). 


Figure  II  - Maillage  libre 


Figure  12  maillage  regie 


Figure  13  maillage  hybride 


Figure  14  Visualisation  des  contraintes  dans  la  bride 

La  figure  suivante  montre  revolution  de  la  contrainte  de  Von  Mises  maximale  atteinte  en  fond  de  trou,  en 
fonction  du  rayon  de  fond  de  feston. 


Contrainte  en  Fond  de  Trou 


Figure  15  Evolution  des  contraintes  maxi  en  fonction  du  rayon  de  feston 
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Ces  courbes  permettent  de  constater  que  les  resultats  obtenus  sont  tres  sensibles  au  type  de  maillage  employe. 

Le  niveau  moyen  de  la  contrainte  est  different  pour  chacune  des  courbes.  Cela  s’explique  par  le  fait  que  les 
mailles  sont  de  tallies  differentes  suivant  le  type  de  maillage.  Les  pentes  sont  aussi  tres  differentes.  Cette 
difference  s’explique  par  la  variation  de  I’element  le  plus  contraint.  Cette  variation  est  d’autant  plus  grande 
que  Ton  utilise  un  maillage  libre. 

Cette  etude  montre  qu’en  I’absence  de  perturbateur  de  maillage  une  optimisation  ne  pent  etre  realisee  sur  un 
maillage  libre.  En  effet,  les  sensibilites  sont  fausses  et  empechent  les  algorithmes  d’ optimisation  de  converger. 

Cependant  la  complexite  des  pieces  etudiees  rend  parfois  delicat  la  realisation  d’un  maillage  regie.  Combinant 
les  avantages  de  I’une  et  de  I’autre,  la  strategie  de  « maillages  hybrides  » est  un  bon  compromis  pour  etre 
appliquee  dans  une  chaine  d’optimisation. 

Conclusion 

Les  deux  exemples  proposes  illustrent  qu’en  quelques  annees  les  techniques  d’optimisation  ont  evolue  de 
fagon  importante  et  permettent  desormais  de  traiter  des  problemes  d’optimisation  de  forme  de  plus  en  plus 
complexes.  Ainsi,  il  est  possible  d’utiliser  tons  les  types  d’analyse  mecanique  et  egalement  de  coupler  les 
logiciels  de  CAO,  les  logiciels  de  maillage  et  les  solveurs  elements  finis. 

Neanmoins,  ces  exemples  illustrent  egalement  certains  dangers  dans  I’utilisation  des  logiciels  sans  reellement 
en  maitriser  tons  les  aspects  (calcul  des  sensibilites  errone  par  une  qualite  de  maillage  insuffisante,  absence  de 
pertubateur  de  maillage,  fonctions  en  « escalier » non  continues  et  conduisant  a des  sensibilites  nulles, 
problemes  des  minima  locaux,  problemes  multi-objectifs,  approximation  insuffisante  des  surfaces  de  reponse, 
etc  ). 

Nous  avons  note  d’autres  freins  d’utilisation  de  I’optimisation  qui  sont  lies,  non  pas  a I’optimisation  en  elle- 
meme,  mais  a la  difficulty  de  realisation  des  modeles  parametres  (CAO  ou  maillage  automatique 
completement  associe  a la  CAO)  et  a la  difficulty  de  poser  correctement  un  probleme  d’optimisation. 

Nous  travaillons  pour  lever  ces  quelques  imperfections  et  il  nous  apparait  dgalement  ndcessaire  de  former  et 
d’aider  les  ingynieurs  d’un  bureau  d’dtudes  afin  qu’ils  acquierent  « un  niveau  de  culture  en  optimisation  » . 

Pour  conclure,  nous  sommes  convaincus  que  la  facility  d’utilisation  des  futurs  logiciels  de  CAO  pararndtrys 
comme  CATIA  v5  et  I’intygration  de  la  plupart  des  logiciels  de  maillage  dans  ces  environnements  vont 
permettre  de  multiplier  de  fagon  importante  les  dtudes  d’optimisation  et  de  rendre  plus  industrielles  les  dtudes 
multi-disciplinaires. 


Paper  #13 

Discusser's  name:  Kirit  Patel 
Author’s  name:  Mr.  Girard 

Q:  In  your  pre-optimization  process,  you  spent  some  time  in  the  design  of  the  disc  you  showed.  How 
much  improvement  in  design  time  has  resulted  from  your  new  optimization  method? 

A:  The  optimization  techniques  allow  us  to  explore  a much  larger  space  in  the  design  area  for  the 
same  amount  of  time  as  traditional  techniques.  Consequently,  the  main  improvement  is  not  in  design 
time,  but  in  quality  of  the  final  design.  This  results  in  direct  reduction  of  time  and  cost  in  the  global 
design  process:  it  helps  to  reach  the  'Best  Design  at  first  time'  and  reduce  non-quality  & re-design, 
therefore  reducing  cost  and  delays  on  a program 


